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Abstract 

For any N > 2 and := (ai,--- ,aAr+i) € (0,oo)^+^, let be the corre¬ 
sponding Dirichlet distribution on A := {x = (xi)i<i<Ar G [0, l]'^ : Yi<i<7V — l}- 
We prove the Poincare inequality 


1 r ^ 

< -- I {(l - X + f ^ 

^ l<i<N n=l 


and show that the constant is sharp. Consequently, the associated diffusion 

process on A converges to in at the exponentially rate aw+i. The 

whole spectrum of the generator is also characterized. Moreover, the sharp Poincare 
inequality is extended to the infinite-dimensional setting, and the spectral gap of the 
corresponding discrete model is derived. 
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1 Introduction 

Let G N. For any a = (ai,-- - ,aAr+i) G (0,oo)^'''^, the Dirichlet distribution with 
parameter a is a probability measure on the set 

:= jx = (a:i)i<i<Ar e [0,1]^ : ^ x* < l| 

l<i<Af 


with the density function 


lll<i<Ar+l l<i<Ar 

where |a:|i := X]i<i<Ar |xi| for a: G Obviously, identihes to the distribution 

fl^^+^\dx,dy) := y^^\dx)6i_\^\^{dy) 

on the space 

y(Af+i) _ G [0,1]^+^ : y+ |x|i = l|. 

The Dirichlet distribution arises naturally in Bayesian inference as conjugate priors for 
categorical distribution and inhnite non-parametric discrete distributions respectively. They 
also arise in population genetics describing the distribution of allelic frequencies (see for 
instance [lliniEn]). In particular, for a population with A+1 allelic types, a:j(l < i < A+1) 
stands for the relative frequency of the Tth allele among A + 1 ones. 

The Dirichlet distribution possesses many nice properties. We will use the following 
partition (or aggregation) property of for a G (0,oo)^’''^. Let (Ai,..., Aat+i) have 

law let Ai,A 2 , ..., Ak+i be a partition of the set {1,2,..., A + 1}, and set 

Yj = ^ ^ Ar, fdj = ^ ^ ftr) j = 1,k 1. 

r&Aj r&Aj 

Then (Ai,..., Yfc+i) has law with parameters [3 := {j3i,..., jdk+i) G (0,cxo)*^+^. We 

would also like to recall the neutral property of the Dirichlet distribution ([!]). For (Ai, • • • , Atv) 
having law we dehne 

Di = Ai, Ui = - -— 2<i<N. 

Then Ui is a beta random variable with parameters {ai, «*+! + ... + aAr+i) and Ui,... ,Un 
are independent. This leads to the following representation of the random variable with law 

A^-1 

(Ai, A2, ..., A^) = (Di, U2{1 - Di),..., [/^ n (1 - • 
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A well known construction of the Dirichlet distribution is through a Polya’ urn scheme 
([!]). More specihcally, consider an urn containing A^ + 1 balls of different colors labelled by 
+ The initial mass of the ^colored ball is ai. Balls are drawn from the urn 
sequentially. The chance of a particular colored ball being selected is proportional to the 
total mass of that colored balls inside a urn. After each selection, the ball is returned with 
an additional ball of same color and mass one. The relative weight of different colored balls 
inside the urn will eventually converge to a Dirichlet vector (Xi,X 2 ,... jX^v+i). 

To simulate the Dirichlet distribution, several diffusion processes with this distribution 
as the stationary distribution have been proposed and studied. The Wright-Fisher diffusion 
(see [El [iHl Uni [22]) is a diffusion approximation to the Wright-Fisher Markov chain model in 
population genetics. The evolution mechanism involves mutation and sampling replacement. 
It is reversible with respect to the Dirichlet distribution. Exploring the property of right 
neutrality, a GEM diffusion is introduced in m and studied further in [12]. This is also 
a reversible diffusion with Dirichlet distribution as the reversible measure. A key problem 
of the study is to estimate the speed for the diffusion process to converge to the Dirichlet 
distribution. 

The inhnite-dimensional generalization of the Dirichlet distribution is Ferguson’s Dirich¬ 
let process ([E]). It is a random atomic probability measure characterized by the property 
that its restriction on any hnite partition of the state space is a Dirichlet distribution. The 
masses of the atoms follow the GEM distribution (see m)- The inhnite-dimensional gener¬ 
alization of the Wright-Fisher diffusion is the well known Fleming-Viot process with parent 
independent mutation ([lll[9]) which has an unlabelled version, the inhnitely-many-neutral- 
alleles model ([8]). The Dirichlet process is the reversible measure of the Fleming-Viot 
process. Various functional inequalities have been studied to investigate the efficiency of 
these processes in approximating the Dirichlet processes (0. IZZl). In particular, it was 
shown in [22] that a logarithmic Sobolev inequality holds for the Fleming-Viot process if and 
only if the state space has hnite dimension. The Wasserstein dihusion studied in [21] and [5] 
is closely related to the Dirichlet process with state space [0,1]. However, the exact conver¬ 
gence rate is not yet known for these processes to approximate the Dirichlet distribution. 

In this paper, we will construct and study a new class of processes with Dirichlet distri¬ 
bution and Dirichlet processes as reversible measures. In comparison with the Wright-Fisher 
Markov chain model where sampling replacement occurs between any pair of individuals, our 
hnite dimensional Markov chain model only allows sampling replacement between individ¬ 
uals in one group and individuals in another group. This reduced sampling scheme bridges 
the gap between the independent systems and the pairwise sampling models. A complete 
understanding of these models will provide a whole picture of the roles played by diherent 
evolutionary forces. 

The main contributions of this paper are the explicit identihcation of the whole spec¬ 
trum of the hnite-dimensional dihusions, the construction of an inhnite-dimensional dihusion 
process with Dirichlet process on a countable state space as the reversible measure, the estab¬ 
lishment of sharp Poincare inequalities for both the hnite and inhnite dimensional dihusions, 
and the construction of Markov chain models. In particular, we found the exact exponential 
rate of the underlying dihusion processes and Markov chains to converge to the Dirichlet 
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distribution and its infinite-dimensional generalization. 

An outline of development of the paper is as follows. In Section 2, we collect the main 
results. The sharp Poincare inequality for the finite-dimensional diffusion is proved in Section 
3. The key step in the proof is to link the eigenvalues of the diffusion generator with a finite 
matrix. The whole spectrum is obtained in Section 4. Section 5 contains the construction of 
the infinite-dimensional diffusion and the establishment of the corresponding sharp Poincare 
inequality. Finally, in Section 6 we introduce the discrete Markov chain model involving 
immigration, emigration and sampling, which approximates the diffusion process solving 
(12.Ij) . This not only provides the genetic link to the diffusion model but also opens the door 
for the study of models with other sampling mechanisms. 


2 Main Results 


The diffusion process studied in this paper first appeared in [121 (2.44)] (see also [2]) and 
solves the following SDE on 


El 


(2.1) dX,{t) 


{a^{l-\X{t)\)-aN+iXi{t)}dt+,/2{l-\X{t)U)Xi{t) d5,(f), 1 < i < N, 


where B{t) := ■ ■ ■ ,B]\f{t)) is the d-dimensional Brownian motion. 

We will show that the Markov semigroup associated to fl2.1l) is symmetric in L^(yU,o^^); 
that is. 


E2 


( 2 . 2 ) 


ja(^) Ja 


holds for 

L[^\x) := ^ {xn{l-\x\i)dl +{an{l-\x\i) - aN+lXn]dr^ 

l<n<N 

being the generator of P“, where and dn := So, (Lq,, (^^(A^^))) is closable in L^(/ii^^) 
and its closure (Lq, ^(Lq,)) is a negative definite self-adjoint operator. 

Moreover, since 


N 

= {fL\^'^g + gL^^'^f){x) + 2{l- |x|i) ^x„{(a„/)(«9„^)}(x), 

n=l 

(12.2p implies the integration by parts formula 


E2 = 


(2,3) 


IaW 


r ^ 

fL[^'>gdn‘'f'> = |(1 - \x\i)'^xn{{dnf){dr,g)}{x)^g.^f\dx) 


n=l 


=:4^\f,g), f,geC\A^^^). 


Therefore, is closable in L^(/ii^^) whose closure is a sym¬ 
metric Dirichlet form on and it is easy to see that this Dirichlet form is associated 

to the Markov semigroup P". 
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Tl.l 


Finally, the spectral gap of is characterized as 

gaplLW) = inf {<rW(/,/) : / e >(/) = 


= 1 


It is known that when = 1 we have gap(Li^^) = ai + 02 , see e.g. [22]- So, in the following 
result we only consider N > 2. 

Theorem 2.1. Let N >2. Then P" is symmetric in and its generator has spectral 

gap gap(LL^^) = oat+i. Consequently, P“ converges to exponentially fast in : 

II p“ _ ;,(-^)|| < p-«fv+n + > n 

IIPt Ta IIl 2(^W) P e , c ^ u, 

and the sharp Poincare inequality for is 

f 6 = 0 . 

Oat+I 

Next, we extend this result to the inhnite-dimensional setting. Consider the inhnite- 
dimensional simplex 


:= |x e [0,1]^ : \x\i = '^Xi < l|, 

i=l 


which is equipped with the L^-metric \x — y\i. Let a G (0, 00 )^ with |q;|i = 
and let Ooo > 0 which refers to OAf+i in the hnite-dimensional case as N ^ 00 . Let 


^ I Oi, ■ ■ ■ , On— 1 ) 


, ^ Oi, Ooo) e (0, oo)’"+\ n>l. 


i>n 


Then for any n > 1, 


•=,dxn) Yl 5o(dxi 


2=n+l 


T1.2 


is a probability measure on We will prove that when n ^ 00 these measures converges 

weakly to a probability measure on A’'°°i, which is the inhnite-dimensional generaliza¬ 

tion of Dirichlet distribution with parameters {a, Ooo). The following result extends Theorem 
l2.1l to the case for N = 00 . 

Theorem 2.2. Let a G (0, cxo)^ with |a|i < 00 and let > 0. 

(1) The sequence {pa”ooo}n>i converges weakly to a probability measure on A(°°b 

(2) The form 

p 00 


n=l 


5 







is closable in whose closure is a symmetric Dirichlet form. The generator 

of the Dirichlet form satisfies C and 

OO 

-^Soo/(^) = “ Mi)9lf{x) + {an(l - |a:|i) - a^Xr,}dnf{x)'^, f E ^C‘^. 

n=l 


(3) The generator has spectral gap gap(La“^^) = aoo- Conseguently, the associated 

Markov semigroup converges to exponentially fast in Lf{pi'^aao) '■ 


\pa,a^ _ Aoo) II ^ -a^t +>n 


and the sharp Poincare ineguality is 


/), / e ^C-*,f,>“>„(/) = 0. 


Finally, the next result shows that the diffusion process generated by is the weak 

limit of the -diffusion process as n ^ oo, where 

n n n 


T1.3 


2 = 1 


2=1 


2=1 


For any x G and T > 0, let be the distribution of the diffusion process generated 
by L^a}aoc with initial point x^'^^ := (xi,-- - , x„_i, Xj). Embedding A^"-^ into A(°°) by 

setting Zj = 0 for 2 ; G A*^”) and i > n -|- 1, we regard as a probability measure on 
Dt := C([0,r]; Aoo) equipped with the uniform norm ||^||i,oo := sup^gioT’] |^(t)|i. 

Theorem 2.3. For any x G and T > 0, converges weakly to a probability measure 

P^ on Dt- Moreover, P^ solves the martingale problem of for any f G the 

coordinate process X{t){u) := u{t) and the natural filtration '■= (xipAs '■ s G [0,t]), 

f(X(t))- f V:Aj(X(s))ds, te[0,T] 

Jo 

is a martingale under P^ ■ 
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Proof of Theorem 


2.1 


We hrst prove fl2.2p which implies the symmetry of P“ in Lfi{yi'a'^). Since smooth functions on 
AW are uniformly approximated by polynomials up to second order derivatives, it suffices 
to consider f,g E A^oo, the set of all polynomials on AW. Let 

= Xn{l - \x\i)dl + {an(l - |x|i) - aN+iXn}dn, l<n<N. 
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E3 


Then fl2.2p follows from 
(3.1) 


/ (n n / (n n 

r^, and simply 


"(N) 


l<i<N 


l<i<N 


l<i<N 


for Pi,qi E Z+, 1 < z < A^. Letting Pn+i = Qn+i = 0 and C = 
denote xn+i = 1 — |a:|i, we have 


l<i<Af 

r(|«|i) 


ni<i<Ar+i J 


/ (n n 

JA(^) i<j<^ 


l<i<N 


= C 


/a(^) 




l<i^n<N +1 


CcjriS i^Qn H“ ^) / ( 1 T 

I Jaw 


l<i^n<N-\-l 


^Pi+Qi+ai 1 2^^ 


ttW+l 


Iaw 


n 


X' 


Pi+qi+oii-1 


l<i<N+l 


dx 


ni<i^n<V+l r(Oi + Pi + Qi) 
l<i<V+l(‘^* + Pi + Qi)) 

X + ttn - l)r(a 7 v+l + l)r(Pn + “ 1 ) “ Ov+ir(ajv+l)r(p„ + + ttn)) 

CT(«7v+i + 1) ni<i^n<v+i r(ai + Pi + ?i) 


r(E 


l<i<Ar+l 


>i +Pi + qi)) 


PnqnJ^ijPn T qn T Q?n 1)) 


where the last step is due to the identity r(s+l) = sr(s) ,s > 0. Since the result is symmetric 
in (pn,gn), it implies (EH]). 


For any d eN, let be the space of all polynomials in whose total degrees are less 
than or equal to d. Let ^o.d = {/ ^ ■ d^a\f) = 0}. It is well known that '■= ^d>i^d 

is dense in so that ^ 0,00 := ^d>i^o,d is dense in 

^0 := {/e = 0} 


under the Sobolev norm ||/||i ,2 := \/ d^a\P) + f) ■ 

To characterize gap(Li^^), we make the spectral decomposition of in terms of the 
degree of polynomials. Obviously, every ^o,d is an invariant space of L^\ Let 
and 

= {/ e ^o,d ■ = 0 for all g E ^d-i}, d>2. 


Then, by the symmetry of in every ^d is an invariant space of as well. 

Thus, letting ^ •^d be the orthogonal projection with respect to the inner product 
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= tt.lW/, d>i,fej^oo 


in 


we have 


PI 


(3.2) 


Therefore, to characterize the spectrum of it suffices to consider that of the 

restriction of on ^i, for every i > 1. 

Let d > 2. To characterize the spectrum of let 


K,= {k = {h,--- ,k^)eZ^: ^A;(z) = d}. 

l<i<V 


For any k G Kd, let = ni<i<v ■ Then 
0 (3.3) ~ ^ (cfc)fceKd e M^A. 

^ k&Ki k&Ka ^ 

We dehne the Kd x iL^-matrix Md by letting 

{ daN+i + Y.i<n<Ni.^n + an-l)kn, H k = k', 

{kn + an)ikn + 1), if /c' = /c + e„ - 6 ^, 1 < ^ 7^ m < N, 

0 , otherwise, 


where {e„}i<„<Ar is the canonical orthonormal basis on . We hrst identify eigenvalues of 
with those of Md- 

I LI I Lemma 3.1. For any d> 2, X is an eigenvalue o /—if and only if it is an eigenvalue 
of Md- Conseguently, —> {daN+i)!^^, where is the identity operator on £2d- 

Proof. (1) Let A be an eigenvalue of on By fl3.3l) and (13.21) . there exists 0 7 ^ c G 

such that 


W1 


(3.4) 


Obviously, 


Ck{L^^^x^ - Tid-iL^^^x'^) = -A Ck{x^ - Tld-ix'^). 
k£Kd k&Kd 


iW^k _ ^ 


l<n<N 


l<n,m<N 


l<n,m<N 


X 


f ^ l)x 






l<n<N 


n,m<N 


l<n,m<N 


^ ^ kni^kji 1 ) 

l<n,m<N 


X 


k-en+en 


+ ankr 


X 


+ dAT+i ^ ^ k 
l<n<N 

k-e„+e^ + . 


l<n,m<N 










By the change of variables k' := k — Cn + e-m, 'we obtain 

kGKd l<n,m<N 

= '^ Ck + X] OinknX^ 

kGKd l<n^m<N k^K^ l<n<N 

^ ^ ^ ^ ^k'+en—em^ni}^ ^ ^ ^k ^ 

k^Kfi l<nj^m<N k^K^i l<n<A^ 

^ ^ ^ ^ ~\~ l)C/^-j-e^—em^ “ 1 “ ^ ^ ^ ^ ■ 

kGKfi l<n^m<N k^K^ l<n<A^ 

Similarly, 

k£Kfi l<n,m<N 

^ ^ ^ ^ kni^kn -|- l)c/^-(_e^— em,^ ^ ^ ^k ^ ^ ^n(^n l)x . 

k^Kfi l<n^m<N kGK^ l<n<A^ 

Combining these together leads to 


W2 


(3.5) y] CkL^^'^x^ = + <^n5n)a:^ 

fceiCd keKd l<n<N 

Substituting this into fl3.4p . we arrive at 


Md{k,k')ck'X^. 

k,k'£Ka 


iMdc)kX^ = ^ + P<i-i(a^) 

fce^d fc€Kd 

for some pd_i G J^d-i- Therefore, M^c = Ac, i.e. A is an eigenvalue of Md. 

(2) On the other hand, if A is an eigenvalue of Md, then there exists c G \ {0} such 
that MdC = Ac. Let 

f{x) = y] CfcX^ - Tid-i y CfcX^. 

fceiCd fcGJCd 

It follows from M^c = Ac and (13.51) that 

L'Pf = Pi-1 - A/ 


holds for some p^-i G l^d-i- Since / G ^d which is orthogonal to ^^d-h this and (13.21) 
implies 

Lf^f = (1 - nd-i)Li^^f = -A(l - nd-,)f = -A/. 

So, A is an eigenvalue of on ^d- 

(3) Finally, since eigenvalues of —L^a'^ are nonnegative, (2) implies that eigenvalues of 
Md '■= Md — daN+ilR^xKa is larger than or equal to —daN+i. On the other hand, from 
the dehnition of Md we see that Md does not depend on av-i-i- So, letting ckat+i J, 0 and 
noting that Md > 0, we conclude that eigenvalues of Md are non-negative. Therefore, 
eigenvalues of Md are larger than or equal to da^^i. Combining this with (1) we obtain 
—> {daN+i)I^^. □ 
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Proof of Theorem \2.1[ By Lemma 13.11 it suffices to prove that the smallest eigenvalue of 
—is aiv+i. To this end, we take 9i = (%)i<j<v ^ M^(l < i < N — 1) such that 


N 


N 


^ ^ 9ikC^k 0) ^ ^ 9i}~9 f — h J — ^ f- 


k=l k=l 

So, {9i}fLi is a basis of Let 

N 


1 < f — 1; 

1=1 
N 

un{x) = - 


fe=i 


a 



a := o 

|o|i 



V 

E' 

fc=i 


We intend to prove that {ni}i<i<Ar is an orthogonal basis of with respect to the inner 
product {f.,g)^a'’ := pT^\fg) = jf g<lg'a\ and = -|a|iWv while = 

—aN+iUi ioT 1 < i < N — 1. Thus, the smallest eigenvalue of —is a^+i- 
It is easy to see that 

(X.) := / x.^r'(dx) = rnTinr/\ = FT’ 

Jaw r(|a|i + l)r(Q!i) |«|i 




r(a;)r(a;j + 2 ) cri(cii + 1 ) 




r(|«|i + 2)r(aj) |a|i(|Q;|i + 1) 

r(Q;)r(«j + l)r(«j + l) 


r(|a|i + 2)r(aj)r(aj) |a|i(|a|i + 1) 


1 < 1 - 1 ; 


Then 


N 


= TW X! 1 < * < ^ - 1 ; 

\n/ \ ^ 


Q; h 


N 


k=l 




ai a 


= 0 . 


i=l 


oh Ct 1 


So, {ui}i<i<N C Moreover, for 1 < z 7^ j < A^ — 1, 

1 




|ct|i(|ti|i + 1 ) 


|ci|i(|<a|i + 1 ) 


^ ^ 9ik9 jk^kiS^k T f) T ^ ^ 9■ik9 jlOLkOLi 

l<k^l<N 


l<k<N 


^ ^ 9ikO!k^ ^ ^ 9jlCXl T ^ ^ 9^i^9ji^CXk r 0, 

1 <;<A'' T<ri.<rM ' 


l<k<N 


l<k<N 


10 
















and for any 1 < i < N — 1, 

ljS^\uiUN) = 'Yl 

l</cj'<W 

'Y + + 


|®|i(|ck|i + 1 ) 


|Q!|l(|ck|l + 1 ) 




^ ^ (^ijOijOlk “ 1 “ 


q;|i(|q;|i + 1 ) 
1 




l^kJ<N 

\N—1 ■„ „ „f TmA'’-! 


Ct|i(|Q!|i + 1) 


''Y^ ~ O' 


i<j<N 


Since is a basis of we have 

dim span{nj :l<z<n — 1} = A^ — 1 = dim 
In conclusion, {ui}i<i<N is an orthogonal basis of 


Finally, we have 


N 


L^^'>Ui{x) = ^(oj-xat+i - aN+iXj)9ij = -aN+iUi, 1 < i < iV - 1, 

i=i 


and 


N N N 

L^^\n{x) = Yji^j^N+l - OlN+lXj) = -\a\iY^^j + Y^^^ ^ -\<x\iUn{x). 

j=l j=l j=l 

Therefore, the proof is finished. 

(TV) 


□ 


4 The whole spectrum of La 

For d G Z+, let be the space of homogeneous polynomials of total degree d in the variables 
Xi, ..., xat. Denote by Hd the natural projection from to which only keeps the d- 
homogeneous part of a polynomial. Let = {TTdLL'^)\j^^ be the restriction of the operator 

to Md and denote —l^d its spectrum, seen as a multi-set (namely with multiplicities). 
From the above considerations, the spectrum A of is equal to yJd&+-^di as a multi-set. 

We can write 

j{N) _ I I j{N) f{N) 

^a,d - -\ • n^a,d - <^N+l^a,d^ 

where : M’d —t d^d-\ and : M'd ^d are respectively the restriction to M’d of the 

operators 

■= Y i^nd'^ + andn), Y ^nOn- 

l<n<A^ l<n<N 

The crucial point of the previous decomposition is that = dljfr^. Denote by the 
spectrum of | • we thus have 

Ad = Arf -|- daN+i- 

Note that Aq = Aq = {0}. The next result enables to compute by iteration Ad for all d G Z+. 
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pro2 


Proposition 4.1. For any d G Z_|_, we have 


Ad+i = {2d+a + Ad) U {0[C{N, d+1)- C{N, rf)]}, 


where {0[/]} is the multiset with 0 repeated I times, for I G Z_|_ (more generally [1] will stand 
for the multiplicity 1), and where C{N,d) is the dimension of J^d, namely 


C{N, d) 


ri"} 


Proof. Consider A G A^+i 


and let ip G J^d be an associated eigenvector (non-zero). We have 


= V- 


Since belongs to Jifd, there are two possibilities: either A = 0, or (p 

some G J^d such that 


*W 


(4.1) 


f{N) 






if; for 


We consider the latter situation, since the former case leads to the multi-set {0[C'(A^, d -|- 
l)-C(iV,d)]}. 

We compute at point x that 


+W2 


^S+l(l • M -11 + 2 

l<n<Ar 

(4-2) = + if ^ 0,1 + 2 ^ Xndnfj 

l<n<A^ l<n<N 

= + (d + ‘2d)fj. 

So, it follows from fl4.ip that A — d — 2d is an eigenvalue of the operator | • namely 

belongs to A^. Thus, 


Arf+i c (2d + d + Ad) U {0[C'(A^, d+1) - C{N, d)]}. 

On the other hand, if A' G A^ then | ■ = X'fj for some 0 7 ^ G Then fl4.2p 

implies 

-bi^j+id • 11-0) = I ■ 11+^0 + (d + 2d)0 = (A' + d + 2d)0. 

Therefore, A' + d + 2d G A^+i; that is, A^+i D (2d + d + A^). Then the proof is hnished. □ 

The previous arguments amount to an iterative construction of the eigenvectors: for any 
d G Z+, let ^d be the set of eigenvectors of and ^d be the kernel of Then we have 


Vd G Z+, <^d+i — ^d+i 0 yN^d- 
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Indeed, in the above proof, functions (p G ^d+i of the form yjq'ip with ip G are associated 
to eigenvalues of the form d + 2(i+ A, where A G Ad- From Lemma [3.11 we know that A > 0, 
so that d + 2(i + A > 0 and ip does not belong to the kernel of L^^i- Conversely, we have 
seen that all the other eigenvectors belong to the kernel of Thus we get the following 

characterization of the kernel of : it consists exactly into the eigenvectors of which 
don’t admit ?/Ar as a factor. 

Note that ^d is also the set of eigenvectors of L^al- To get the eigenvectors of our initial 
operator L, we construct by iteration on d G Z+ the following subsets ^d of First we 
take ^0 := ^0 = Next, if has been constructed, then for any / G ^d+i, there exists 
a unique Qf E S^d such that / + (?/ is orthogonal to S^d iu A^(;u). Then we define 

^d+\ '■= {f A gf j E ^d+i}- 

The set of eigenvectors of L is U^GZ+^d- 

From Proposition 14.11 it is possible to parametrize the spectrum A of —L in the following 
way. Let be the set of elements of the form (/ci, A: 2 ,...,/cr, ^r+i), where r G Z+ and 
0 < /ci < /c 2 < • • • < < K+i- Define a mapping K : PXf —)■ A via 

'ik := (fci, A: 2 , ..., kr, K+i) E J^, K{k) := 2{ki + ■ ■ ■ kr) + ra + K+ia^+i- 

Then K is surjective. It is truly one-to-one, if and only if 1, d and oat+i are independent 
when M is seen as a vector space over Q. Let us call this situation generical over the choice 
of the parameters a := (a„)i<„<Ar+i. 

The multiplicities can also be recovered. Consider the mapping D : ^ —)■ N defined by 
D{k)-.= Y,C{NM)A Y. {C'(iV,/ + l)-C'(iV,/)} 

l</<r ,/i;2 )^r} 

for k := {ki, /c 2 ,..., K, K+i) G J^. Then the multiplicity of an eigenvalue A G A is given by 

fcGK-l(A) 

In particularly, generically, we have A = {K{k)[D{k)] : k E J^}. 


5 Proofs of Theorems 12.21 and 12.3 


To prove the first assertion, let W be the L^-Wasserstein distance induced by p(x, y) := \x — 
l/|i on the set of all probability measures on That is, for any E 

W{y,i'):= / \x - y\7i{dx,dy), 

JTTe{ii,u) 7a(°°)xA(°°) 

where ^(/i, zz) is the set of all couplings for y and z/; i.e. tt G (p., z^) if and only if it is a 
probability measure on x such that 

7r(da; x A^°°^) = p(da;), 7r(A*'°°^ x dz/) = z^(dz/). 
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It is well known that the metric W is complete and induces the weak topology on 

see e.g. [3l Theorems 5.4 and 5.6]. So, for the proof of Theorem 12.21 we only need to show 

that {pa”ooo}n>i is IT-Cauchy sequence. 

Proof of Theorem 12.21 (1) To prove that }n>i is a IT-Cauchy sequence, we use the par¬ 

tition property of the Dirichlet distribution mentioned in Section 1. For any n > m > 1, let 
(Xi, • • • , Xn+i) have law By the partition property, (Xi, • • • , X^-i, ^n+i) 

has law So, (Xi,--- , Xm-i,'Ei^^Xi) has law while (Xi,-- - ,X„) has law 

Thus, the laws of (Xi,-- - ,X^_i,Jf;^^^Xi,0,0,--- ,0) and (Xi, • • • ,X„,0,0,-- - , 0) 

are and respectively. Then, by the definition of W and noting that |a|i < oc, 

we have 


limsup sup < 2hmsup sup ^ E|Xi 


m-^OQ n>m+l 


m—>-oo n>m+l 


i=m+l 


2aj 

= hm sup sup 2^ - ^oo - 

m^oo n>m+l ®oo + Z^i=n +1 


= 0 . 


Therefore, n>i is a hF-Cauchy sequence and the proof of the first assertion is finished. 

(2) It suffices to prove 


^ (5.1) 


(/■ 9) = - / (/4“l9) . /. 9 e 

J A(°°) 


For any f,g E ^(7^, there exist m G N and fm,gm E such that 

f{x) = fm{Xi,---,Xm), g{x) = gm{xi,- ■ ■ ,Xm), X G 
So, by the definition of g^ajxoo and using fl2.3p . we have 


^ (5.2) 


f A(°°) 


f l<i<n i=l ' 


n) 

Oio 


Since —)■ weakly, and it is easy to see that 

hm sup IfL^^l^g - fL^'^J^glix) = 0, 


' xeA(°°) 


hm sup 

n->-oo 3.g^(oo) 


Xi(dif)(dig) - (l - Xi(dif)(d,g) 


l<i<n 2=1 


2=1 2=1 


= 0 , 


by letting n —)■ cxo in fl5.2p we prove fl5.ip . 

(3) Finally, as was shown in (2) that the desired Poincare inequality follows by applying 
Theorem 12.II to on then letting n oo. So, gap(Li')^^<^) > a^o- On the other hand, 
let 

u(x) = a 2 Xi — aiX 2 , x E X^°°\ 
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We have 

= {“i(l “ l^li) ~ aoc,Xi}a2 - {a2(l - k|i) - aooX2]cii = -aoou{x), x G 


This implies gap(Li'^^<^) < cioo- In conclusion, we have gap(La'^''^) = a, 


(oo) 


□ 


Proof of Theorem \2.S\. (a) For the first assertion, we only need to prove that {-Pj^}n>i is a 
Cauchy sequence with respect to the L^-Wasserstein distance 

To this end, for any n > m > 2, we construct a coupling of and P^f^ as follows. 

Firstly, let (X-^"^(f))i<j<„ solve the following SDE with 


df 


FY 


(5.3) 


FY2 


dx}^\t) =[«,(! - |XW(f)|i) 

+ ^2(l-lX(-)(t)li)x[^\t) dBi(t), l<z<n-l; 

OO 

dX<r>(t) =[l]aj(l - |V”>(()|i) - a„X>">(()]dt 

j=n 

+ ^ 2 {i-\x(-){t)U)xi^\t) d5„(f), fe[o,r], 

where {Bi{t)i<i<n are independent one-dimensional Brownian motions. Then is the 
distribution of {X^'^\t))t(z[o^T]- 
Next, let 

(5.4) Xl^\t) = Xl^\t) for 1 < z < m - 1, and X^j\t) = ^ ^ [0,T]. 


j=m'‘ 


Then X("^i(0) = and by 

dXf )(f) =L(1 - |X(-)(f)|i) - a^X^^\t) 


df 


+ 1 / 2 ( 1 - |a;(™)(f)|i)xf”*^(t) dPi(f), 1 < z < m - 1; 

OO 

dA'r>W =[^0,(1 - IxWli)!,) - aeoAW(t)]di 


]=m 


+ ^J 2 {l-\x^-^){t)\^)xt\t)dB^{t), fe [o,r] 


where dBmif) := 


■/ 


X. 


(m) 


(t) 


Yl^=m \/dBi[t) is a one-dimensional Brownian motion in¬ 


dependent of (Pj(f))i<j<m-i. Therefore, (X^'"i(f))fe[o,T] has law P, 


(m) 
x.T • 
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Now, by fl5.4p and the definition of Wr, we have 


FY3 


(5.5) 


n 


WriP^^ < E sup = 

t€[0,T] 


E sup Xj^\t). 

j=rn+l 


Let Z{t) = X]j=m+i By (|5.3p we have 

CO n I 

dz(t)<( Y1 E v»(i-i'^“wii)N”’WdB.w. 

j=m+l j=m+l 

So, 


Z{t) < (xj + taj) + Y 

j=l-\-m j=m-\-l ' 


s(l - |XW(s)|i)Xf)(s)di?,(5) =: z{t), t e [o,r]. 


Since Z{t) > 0, Z{t) is a nonnegative submartingale. Then by Kolmogorov’s inequality. 


P 


sup z{t) > )^ <f( sup z(t) > a') < —EZ(r) 

tg[o,T] ' '^tG[o,r] ^ A 


3 ^ Y + A > 0. 

j=m+l 


Since Z{t) < 1, this implies 


1 

E sup Z{t) < A + pf sup Z{t) > )^ < A + — 'Y, {^j + A > 0. 

i6[0T] ^te[o,r] ^ ^ 


Taking A 




{xj + CKjT), and combining with (15.dp . we obtain 


lim sup WriP^'^T^ 

m—)-cx) 72>m+l 



CO 


< 2 lim 

m—^oo 




jf=m+l 


0 . 


Therefore, the first assertion is proved. 

(b) Let / e We have f{x) = f{xi, • • • , Xm) for some m > 1 and / G 

For the coordinate process X{t), define 

= f(X(t)) - [‘L‘-:lj(X(a))da. n>m,te [O.T]. 

Jo 

Then (Mi”^)tg(o,r] is a Pj ^-martingale; that is, for any 0 < s < f < T, and any bounded 
Lipschitz continuous function g on Qt measurable with respect to ^s, 


WF4 


(5.6) 


M<">(i)(c..)g(c..)dFi5= / M<")(g)(a,)g(i.,)dPiy 
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We intend to prove the same equality for ^tid 

Af<“>(i) - nX(t)) - f‘ Li~JJ(X(s))d.s, t 6 |0,T]. 

Jo 

By an approximation argument, we may and do assume that / G In this case, 

is bounded and Lipschitz on VLt uniformly in n > m and t G [0,r]. Since g is 
bounded and Lipschitz on VIt as well, there exists a constant C* > 0 such that 

|(Af(">(i)9)(0 - (Af<">g)(()(q)| < C||4 - 1)111,00, n > »,«,>)£ (Jt.* £ |0.r]. 
Therefore, 




< CWt{p!,% n>m,te [0,T]. 


Combining this with flS.bp . lim„_yoo 1TT(-Pi T)t^) = 0) hm„_yoo and noting 

that g}„,>m are uniformly bounded, we conclude that 


' fix’ 

= lim 

n—>-oo 




' Qx 




(oo) 


< 2 ^ lim sup WT(PiS, ) = 0 - 


Then the proof is finished. 


□ 


6 A Discrete Model 

For any iV>l,M>iV + l, consider a population of M individuals of + 1 different 
types. Divide the population into two groups: group I of types 1,..., iV and group II of type 
+ 1. Focusing on group I and treat group II as outsiders or external sources. Initially 
the number of type i individuals is m*, i = 1,..., A^ + 1. The group I evolves as follows: a 
type i individual independent of all others will wait for an exponential time at rate ccat+i 
and at the end of the waiting emigrates to the outside becoming type A^ + 1 ; an outsider will 
independently wait an exponential time with rate and immigrate to group I becoming 
type i] in addition to emigration and immigration, each couple between a type I and a type 
II waits for an exponential time with rate 2 and when the clock rings, either the group I 
individual moves out becoming an outsider or the group II individual moves in becoming 
the type of the selected individual in group I. 

Let X{t) = ..., Mjsiit)) denote the relative frequencies of individuals of dif¬ 
ferent types in group / among the whole population at time t. For a G (0, we 
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construct X{t) as a multivariate Markov chain with generator 



N 


mY, 




2=1 


M 


f{^- 77 ) “ /(^) + - l^li) /(^ + 77 ) - /(^) 


M 


N 

+M^ ^(1 - \x\i)xiy(^x -^)+f{^ + ^)- 2/(a;)}, / e 
2 = 1 


for X G := |x G : |^c|i = — l}; where e, is the unit vector in the ith 

direction. Letting M —)■ cxo and x —)■ y G one gets ^iS/(x) ^ L^^^fiy). 

We will see that the hnite Markov chain generated by £^^1 ^ reversible with 

respect to the probability measure 


W (^\ [(^N+i]m{1-\x\i) TT [oii]Mxi 

) ■ Z{M{1 - |x|i)}! (Mxi)!’ 

where [a]m '■= Y\a=o^^ + z) for a > 0 and m > 1, [ajo := 1, 


X G A 


(TV) 

M 


and 


Z ; = 


[Q;Ar+l]M(l-|a;|i) TT 

{M(l-|x|i)}! 


is the normalization. Moreover, for N >2, has the same spectral gap as 

Theorem 6.1. Let N >2. The Markov chain generated by is irreducible and reversible 
with respect to Moreover, £^^1, spectral gap ow+i in 

Proof, (a) Denote 7i = H for 1 < z < A. For any x, z/ G A^\ let 

r UxiUN+i + M‘^Xi{l - |x|i), if z/ = X - 7i, 1 < z < A; 

qx,y = < aiM{l - |x|i) + M‘^Xi{l - |x|i), if z/ = X + 7i, 1 < z < A; 

1^0, otherwise. 

We have 

(lxy{f{y) - fix)}, X G A^^y\ 

Since q^^y > 0 when x,y E Aj^J^ with z/ = x ± 7* for 1 < z < A, and Aj^J^ is connected by 
the edges x ^ x ± 7*, we see that the Markov chain is irreducible. 

Next, it is well known that is symmetric in if and only if 


Z2 


( 6 . 1 ) 



x,y E A 


(N) 

M ■ 
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To verify this condition, we only need to consider the following two situations. 

(al) y = X + '^i for some 1 < z < A^. In this case we have M|a:|i < M — 1, and by the 
dehnition of 

h!\yL(y) _ M(1 - \x\i){ai + Mxj) _ ^ 

{aN+i +M{1-\x\i)-l){Mxi + l) Qy^' 

(a2) y = X — for some 1 < z < A^. In this case we have Mxi > 1, and by the dehnition 

r (N) 

of hM,a> 

_ jaN+l + M{1 - \x\i))MXi _ 

{M{l-\x\i) + l){Mxi-l + ai) Qy^' 

In conclusion, fib.ip holds and thus, is symmetric in 

(b) For any d G Z+, consider again the space of all polynomials (in N variables) 
whose total degree is less than or equal to d. For any / G and 1 < z < A^, x 
f{x — 7j) — /(x) and x i-T- /(x + ■ji) — /(x) are polynomials belonging to while 

X I—)■ /(x — 7j) + /(x + 7j) — 2/(x) is a polynomial belonging to ^d- 2 - From the dehnition 
of it follows that ^d is preserved by As in Section 2, we consider for d G Z+, 

(with the convention Note that for d large enough, ^d = {0}, nevertheless, we 

still have 


dez+ 

and the J2d are orthogonal. Furthermore by symmetry of in L^(/z^^), each of the 

^d is preserved by Thus it is sufficient to study the spectral decompositions of the 

restrictions of fo ■^d- But this is exactly the same analysis as in Section 2, because 
there we only used the highest monomials. Indeed, note that for all / G ^d and 1 < z < A^, 

X ^ /(x - 7i) - /(x) + X ^ /(x + 7i) - /(x) - 

are polynomials belonging to ^d- 2 i and 

X ^ /(x - 7i) + /(x + 7i) - 2/(x) - 

belong to ^d-s, where we set = {0} if /c < 0. Thus, for any polynomial / G ^d, the 
polynomials and f have the same highest order term (i.e. the term of degree d), 

so that these two operators have the same spectral gap. 

□ 
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